Deep sub-threshold S production in Ar+KCl reactions at 1J6A GeV 
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We report first results on a deep sub-threshold production of the doubly strange hyperon E~ in a heavy- 
ion reaction. At a beam energy of 1.16A GeV the reaction Ar+KCl was studied with the High Acceptance 
Di-Electron Spectrometer (HADES) at SIS18/GSI. A high-statistics and high-purity A sample was collected, 
allowing for the investigation of the decay channel E~ — > A7r~. The deduced H~/(A + E°) production ratio 
of (5.6 ± 1.2 t\\r) ' 10~ 3 is significantly larger than available model predictions. 



PACS numbers: 25.75.-q, 25.75.Dw 

The doubly strange S~ baryon (also known as cascade 
particle) has, in vacuum, a mass of 1321.3 MeV and de- 
cays (ct = 4.91 cm) almost exclusively into the A-ir~ fi- 
nal state [1]. In elementary nucleon-nucleon (NN) collisions 
near threshold it must be co-produced with two kaons ensur- 
ing strangeness conservation. This requires a minimum beam 
energy of E t hr — 3.74 GeV (^s t hr = 3.25 GeV). In heavy- 
ion collisions, the E" yield was measured at various beam 
energies covered by the RHIC [2], SPS [3, 4] and AGS [5] 
accelerators. Though cooperative processes in heavy-ion re- 
actions allow for particle production below NN threshold, no 



sub-threshold S~ production was observed so far. Predic- 
tions of sub-threshold cascade production at energies avail- 
able with the heavy-ion synchrotron SIS 18 at GSI, Darm- 
stadt, were presented within a relativistic transport model 
[6]. The cross sections of the strangeness exchange reac- 
tions KY — ► 7rS (Y = A, E), which are essential for 5 cre- 
ation below the nucleon-nucleon threshold, were taken from a 
coupled-channel approach based on a flavor SU(3)-invariant 
hadronic Lagrangian [7]. The S~/A ratio was found to 
amount to a few times 10~ 4 , varying with system size and 
beam energy, however, being fairly independent on centrality. 
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At SIS 18 energies, also K production, being a prerequisite of 
the above strangeness exchange reactions, proceeds below the 
NN threshold (E t hr — 2.5 GeV). A similar strangeness ex- 
change reaction like that relevant for S production is found to 
be dominant in sub-threshold K production in heavy-ion colli- 
sions [8], i.e. ttY — > KB (B = N, A). Thus, medium effects 
on strange meson properties, like effective anti-kaon masses 
and hence reduced production thresholds, could strongly in- 
fluence the S yield. However, the authors of [6] found the S 
yield to be more sensitive to the magnitude of the cross sec- 
tions of strangeness-exchange reactions than to the medium 
effects due to modified kaon properties. Generally, the yield 
of multi-strange particles, measured below their production 
threshold in NN collisions, is expected to be sensitive to the 
equation of state (EoS) of nuclear matter. In heavy-ion reac- 
tions, the necessary energy for the production of these par- 
ticles is accumulated via multiple collisions involving nucle- 
ons, produced particles and short-living resonances. The cor- 
responding number of such collisions increases with the den- 
sity inside the reaction zone the maximum of which in turn 
depends on the stiffness of the EoS. 

In this Letter we report on the first observation of sub- 
threshold S~ production in heavy-ion collisions. The exper- 
iment was performed with the High Acceptance Di-Electron 
Spectrometer (HADES) at SIS 18 [9]. HADES, primarily de- 
signed to measure di-electrons [10], offers excellent hadron 
identification capabilities, too [11-13]. 

A 40 Ar beam of about 10 6 particles/s with kinetic energy 
of 1.756A GeV (y/s NN — 2.61 GeV) was incident on a four- 
fold segmented nat KCl target with a total thickness of 5 mm 
corresponding to 3.3 % interaction probability. The beam en- 
ergy is known with a precision of about 1CP 3 . The energy loss 
of the beam particles in the target is estimated to be less than 
0.5 AMeV per target slice. The position resolution of the pri- 
mary (reaction) vertex amounts to 0.3 mm in both transverse 
directions while in beam direction it amounts to 1 .5 mm as ex- 
pected from the finite thickness of the target slices. The data 
readout was started by a first-level trigger (LVL1) decision, re- 
quiring a minimum charged-particle multiplicity > 16 in the 
time-of-flight detectors. The integrated cross section selected 
by this trigger comprises approximately the most central 35 % 
of the total reaction cross section. About 700 million LVL1 
events were processed for the present S~ investigation. 

In the present analysis we identified the A hyperons through 
their decay A — > pn~ . Note that the reconstructed A yield 
includes the decay A's of the (slightly heavier) S° baryon 
decaying exclusively into A and a photon [1] which cannot 
be detected with HADES. Hence, the A yield has to be un- 
derstood as that of (A + S°) throughout the paper. To al- 
low for A selection various cuts on single-particle and two- 
particle quantities were applied. The most important ones 
act on geometrical distances, i.e. i) minimum values of the 
p, 7r~ track distances to the primary vertex (p-VecToPrimVer, 
7Ti-VecToPrimVer), ii) an upper threshold of the p-ir~ min- 
imum track distance (p-7Ti-MinVecDist), and iii) a mini- 
mum value of the A vertex distance to the primary vertex 
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FIG. 1: The p-ir invariant mass distribution. Hatched histogram: 
Scaled combinatorial background produced via event mixing. 



With these conditions we first analysed the invariant-mass 
distribution of proton- n~ pairs (Fig. 1). A clear A signal could 
be separated from the combinatorial background (bg) as de- 
termined via the event mixing technique. The background 
normalization was performed over the given invariant-mass 
range, except a ±15 MeV interval around the A peak. For a 
±3<T mass cut around the A peak, the signal-to-background 
ratio and the significance, defined as signal/ -^/signal + bg, 
amount to 4.6 and 195, respectively. The total A yield for the 
given cuts is N\ = 46, 500 ±600. Fitting a Gaussian function 
to the signal, we obtain a mean value of (1114.2 ± 0.1) MeV 
and a width (a) of (2.7 ± 0.1) MeV. 

Taking this high-statistics A sample, we started the S~ in- 
vestigation by combining - for each event containing a A can- 
didate - the A with those n~ mesons not already contribut- 
ing to the A. The result was a structureless A-ir~ invariant 
mass distribution. Hence, additional conditions were neces- 
sary: iv) a lower limit on the 2nd ir~ (potential S~ daugh- 
ter) track distance to primary vertex (7r 2 -VecToPrimVer), v) 
an upper limit of the distance of the S~ pointing vector w.r.t. 
the primary vertex (S-VecToPrimVer), vi) a maximum value 
of the minimum track distance of the A and the 2nd ir~ 
(7r 2 -A-MinVecDist), vii) a minimum value of the distance of 
the S~ vertex relative to the primary one (S-VerToPrimVer), 
and viii) a window of ±7 MeV around the A mass peak of the 
p-7r~ invariant mass distribution. 

The conditions on the geometrical quantities are summa- 
rized in Fig. 2, where the optimum cut values are indicated by 
arrows. We studied the stability of the signal if more strin- 
gent conditions were chosen. Because of the limited statis- 
tics all other cuts were kept fixed at the optimum values when 
varying a single cut quantity. The dependences on the various 
geometrical distances of experimental data and GEANT [14] 
simulations (see below) are found in good agreement. 

Figure 3 shows the invariant mass distribution of A-ir~ 
pairs after applying all conditions. Indeed, a narrow signal 
shows up on top of a smooth distribution. For an invariant- 
mass window of ±10 MeV (4 bins) around the peak center, 
we find 7V H - = 141 ± 31 ± 25 entries to be attributed to S~ 
with the given statistical and systematic errors. The signal-to- 
background ratio and the significance amount to 0.17 and 4.6, 
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FIG. 2: Relative E~ yields as a function of the cut value of various 
A and E~ geometrical distances (see text, units are mm). The full 
(open) dots display the experimental (simulation) data. The vertical 
and horizontal arrows indicate the chosen cut values and the region 
of accepted distances, respectively. 
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FIG. 3: Top: The same as Fig. 1, but for A-ir~ pairs. Bottom: The 
invariant mass distribution after background subtraction. The full 
line represents a Gaussian fit to the H~ signal. 



respectively. The given systematic error of the signal is due 
to the signal variation for various histogram binnings, back- 
ground normalization regions and mass windows assigned to 
the signal. These systematic variations are also reflected in 
the significance of the signal of about 4 — 6. The full line in 
the bottom panel of Fig. 3 represents a Gaussian fit to the sig- 
nal. The mean value of (1320 ± 1) MeV is well in agreement 
with the PDG value of 1321.3 MeV [1]. Taking into account 
the bias due to the rather large bin size of 5 MeV to be used 
for statistical reasons, the peak width (a) of (4 ± 1) MeV is in 
fair agreement with GEANT simulations which predict for A 
and S~ baryons almost equal values of about 2.5 MeV. In or- 
der to ensure that not a fake signal is selected, we performed 
a A-side-band analysis. No signal was found when choos- 
ing instead of condition viii) a window in the p-ir~ invariant 
mass of 10 < |M p7r - - (M A )| < 25 MeV. Furthermore, when 
dividing randomly the data sample into two sub-samples, the 
S~ sub-yields were found - within errors - compatible with 
half of the above quoted total yield. 

Corrections for the finite acceptances and reconstruction 
efficiencies were deduced from simulations. Thermal A's 
(S~'s), characterized by the temperature parameter T\ (T H -) 
were generated with the event generator Pluto [15]. The ex- 
perimental A rapidity distribution is found slightly broader 
than the thermal model distribution [11, 16]. Consequently, 
in Pluto we allowed also for anisotropic, i.e. longitudinally 
elongated, phase-space distributions. For this purpose, an ad- 
ditional width parameter of Gaussian rapidity distributions, 
a y , is taken into account. The A parameters are chosen such 
that the simulation reproduces both, the experimental val- 
ues of the effective inverse slope parameter at mid-rapidity, 
T e ff,A = 95 MeV, and the rapidity width, a Vi \ = 0.42 



[11, 12]. Since the phase-space distribution of the S~ is 
not known, we investigated its geometrical acceptance for a 
broad range of the transverse and longitudinal shape param- 
eters, i.e. the inverse slope and rapidity width, T" e //. h- = 
(95 ± 25) MeV and a y , E - = 0.34 ± 0.09, respectively. 
Here, the A inverse slope serves as a reference value. The 
lower limit of T e //,H- matches the measured inverse slope 
parameter of K~ mesons [12, 13] being essential for pro- 
ducing S hyperons via strangeness-exchange processes (see 
above), while the upper limit is set by a similar interval above 
the A slope. We assumed the rapidity width of S~'s to be 
larger than the width of thermal S~'s with a temperature of 
95 MeV but smaller than that of A's. This choice is substan- 
tiated by two facts: Firstly, for a thermal rapidity distribution, 
the width approximately scales with the square root of mass 
(cr y = y/T '/m c 2 ), and secondly, the cascade particle may 
carry less longitudinal momentum than the A hyperon, since 
it contains only one light quark which arises from projectile 
nucleons. With the given parameters and their ranges, we cal- 
culated the HADES acceptance (including the branching ratio 
for the decay A — > p7T~) for the A to e acCiA = 0.160 ± 0.009, 
and for the S" to e acc>H - = (9.6±2.i) ■ !0 -2 - The simula- 
tion data are processed through GEANT modeling the detec- 
tor response. The GEANT data were embedded into real ex- 
perimental data and processed through the full analysis chain. 
Relating the outputs after cuts to the corresponding inputs, 
the A and S~ reconstruction efficiencies were estimated to 
e eff;A = (6.1 ± 0.3) • 10~ 2 and e effiH - = (5.5 ± 0.5) • 10~ 2 , 
respectively. We proved our acceptance and efficiency correc- 
tions by extracting the lambda yield [11] which is found to be 
in agreement with existing data [16]. With the above correc- 
tion factors, the ratio of S~ and A production yields can be 
determined. Such a ratio, when derived from the same data 
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analysis, has the advantage that systematic errors cancel to a 
large extent. The ratio is calculated as 



N~- <=acc,A e e ff,A 



A+£° 



N 



A e acc.H- e cff.H- 



(5.6±1.2tJ;?)10- 3 , 



(1) 



where statistical and systematic errors are given, resulting 
from adding the individual ones quadratically. The statistical 
error in (1) is dominated by the 20% error of the S~ signal 
while the systematic error is governed by the stability of the 
signal against cut and background variation and by the range 
of the parameters T s - and a y ^ 3 - entering the simulation. 
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FIG. 4: The yield ratio H _ /(A + E°) as a function of ^/sjviv or 
^/sjvat — ^/sthr (inset). The arrow gives the threshold in free NN 
collisions. The open star, triangles and square represent data for cen- 
tral Au+Au and Pb+Pb collisions measured at RHIC [2], SPS [3, 4] 
and AGS [5], respectively. The filled circle shows the present ratio 
(1) for Ar+KCl reactions at U6A GeV (statistical error within ticks, 
systematic error as bar). Full line: Statistical model for Au+Au [17]. 



when taking - instead of the calculations along the schematic 
freeze-out curve for Au+Au [17] - the optimum input pa- 
rameters (i.e. temperature, charge and baryon chemical po- 
tentials, strangeness correlation volume) following from the 
best fit to all HADES particle yields in Ar+KCl at 1 .76 A GeV 
[11, 13]. Finally, recent predictions within a transport ap- 
proach [6] (soft EoS with incompressibility Kq = 194 MeV 
at normal nuclear matter density) yield a S~/A ratio of a few 
times 10~ 4 , comparable to the statistical model. 

In summary, we observed the production of the doubly 
strange cascade hyperon S~ in collisions of Ar+KCl at 
1.76 A GeV with a significance of about five. For the first 
time, using the HADES detector at SIS18/GSI, this hyperon 
was measured below the threshold in free nucleon-nucleon 
collisions, i.e. at ^/s^vjv — \Jsthr = —640 MeV. Comparing 
the experimental S _ /(A + S°) ratio to the predictions of a 
statistical model and a transport approach, both ones underes- 
timate the experimental ratio. We conclude that I) the condi- 
tions for the applicability of present statistical models might 
be not fulfilled for such rare-particle production in small sys- 
tems far below threshold, and that II) in transport approaches a 
better understanding is necessary of the strangeness-exchange 
reactions conjectured as the dominant process for cascade pro- 
duction below and close to threshold. Other explanations for 
the unexpectedly high S~ yield, like modifications of strange 
hadrons in the nuclear medium, should be investigated. 
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The deduced S~/A ratio (1) may be compared with the 
corresponding ratios at higher energies [2-5]. Figure 4 shows 
a compilation of S~/A ratios as a function of ^snn- The 
displayed data represent the most central 5-10 % of collisions 
of Au+Au or Pb+Pb. At RHIC and SPS energies hardly any 
centrality dependence of the S~/A ratio was observed [2, 3]. 
So far, the lowest energy at which a S~/A ratio is available is 
^~snn — 3.84 GeV, i.e. an excess energy of +600 MeV above 
the NN threshold [5]. The corresponding ratio, measured at 
the AGS at a beam energy of 6^4 GeV, is found to increase 
slightly with centrality. For central (semi-central) collisions 
it is about three (two) times larger than our value. Indeed, a 
steep decline of the S~/A production ratio is expected below 
threshold, where now the first data point is available. This 
allows for comparisons to model calculations. 

The S~/A ratio has been estimated within a statistical ap- 
proach [17]. While RHIC [2], SPS [3, 4] and AGS [5] data are 
well described, the present experimental S _ /A ratio is under- 
estimated (cf. Fig. 4). A similar small ratio is derived [18] 
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